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The harmonic model of atomic nuclear motions is usually enough for multipole 
modelling of high-resolution X-ray diffraction data; however, in some molecular 
crystals, such as l-(2 / -aminophenyl)-2-methyl-4-nitro-l//-imidazole [Paul, 
Kubicki, Jelsch et al. (2011). Acta Cryst. B67, 365-378], it may not be sufficient 
for a correct description of the charge-density distribution. Multipole refinement 
using harmonic atom vibrations does not lead to the best electron density model 
in this case and the so-called 'shashlik-like' pattern of positive and negative 
residual electron density peaks is observed in the vicinity of some atoms. This 
slight disorder, which cannot be modelled by split atoms, was solved using third- 
order anharmonic nuclear motion (ANM) parameters. Multipole refinement of 
the experimental high-resolution X-ray diffraction data of l-(2 / -aminophenyl)- 
2-methyl-4-nitro-l//- imidazole at three different temperatures (10, 35 and 70 K) 
and a series of powder diffraction experiments (20 < T < 300 K) were 
performed to relate this anharmonicity observed for several light atoms (N 
atoms of amino and nitro groups, and O atoms of nitro groups) to an isomorphic 
phase transition reflected by a change in the b cell parameter around 65 K. The 
observed disorder may result from the coexistence of domains of two phases 
over a large temperature range, as shown by low-temperature powder 
diffraction. 



1. Introduction 

When using accurate ultra-high-resolution X-ray diffraction 
data, the most commonly used harmonic model of the atomic 
nuclear motions may not be sufficient for some molecular 
crystals, even for lighter atoms. Therefore, multipole refine- 




ment without modelling anharmonic nuclear motions (ANMs) 
does not lead to the best electron density (ED) model, as 
revealed by peaks and holes in residual maps; these peaks 
arranged in a 'shashlik-like' pattern in the vicinity of the 
anharmonic atoms are an indicator of third-order ANMs 
(Herbst-Irmer et al, 2013; Meindl et al, 2010) and can be 
modelled by introducing Gram-Charlier or cumulant expan- 
sions (Johnson & Levy, 1974). 

Despite the fact that ANMs have been previously discussed 



in the literature (e.g. Kuhs, 1988, 1992), their reliable separa- 

3tion from the static charge-density distribution parameters, 
ACCESS disorder or librations was questioned (Mallinson et al., 1988; 
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Restori & Schwarzenbach, 1996). Although Iversen et al. 
(1999) distinguished anharmonic nuclear motions from static 
electron density features in a thorium complex structure using 
extremely high-resolution (1.7 A -1 ) data from two very low- 
temperature experiments (at 9 and 27 K), Henn et al (2010) 
were able to separate both contributions for lighter atoms 
(namely P atoms) at lower resolution (1.15 A -1 ) at 100 K. 
Birkedal et al (2004) successfully refined the multipolar 
electron density of urea, while Scheins et al (2010) showed 
that ANMs are necessary for the correct description of the 
charge density of a Zn atom. Finally, Zhurov et al (2011) 
showed that neglecting ANMs in the case of hexahydro-1,3,5- 
trinitro-l,3,5-triazine (RDX) results in unrealistic charge- 
density deformation and Laplacian maps in the region of the 
nitro group. For a similar compound, l,3,5,7-tetranitro-l,3,5,7- 
tetraazacyclooctane (HMX), which has a slightly more 
compact crystal structure, the refined ANM parameters were 
statistically significant, however, their effect on the resulting 
charge-density deformation and Laplacian maps was rather 
negligible. 

The effects related to ANMs are visible only at high-reso- 
lution data and the values representing the corresponding 
refined Gram-Charlier coefficients are often hardly statisti- 
cally significant. Correlatively, the agreement factors do not 
improve noticeably upon the introduction of these new 
parameters. Nevertheless, such a physical model considerably 
reduces residual peak heights (Paul, Kubicki, Jelsch et al, 
2011; see Figs. 4 and 5 therein). To avoid possible correlations 
between ANMs and the remaining ED parameters, the former 
ones should be refined first against high-resolution data and 
then by a joint refinement of both anharmonic and electron 
density parameters in the subsequent refinement steps 
(Mallinson et al, 1988). 

Standard resolution crystal structures of numerous 4- 
nitroimidazole derivatives have been investigated in our 
laboratories, with special attention paid to the weak inter- 
molecular interactions present in these molecular crystals 
(Kubicki et al, 2001; Kubicki, 2004a,b; Kubicki & Wagner, 
2007, 2008; Wagner et al, 2007; Wagner & Kubicki, 2007). 
Further investigations of the high-resolution diffraction data 
using the Hansen-Coppens model (Hansen & Coppens, 1978) 
and quantum theory of atoms in molecules (QTAIM; Bader, 
1994) topological analysis were performed for l-phenyl-4- 
nitroimidazole (Kubicki et al, 2002), l-p'-aminophenyl)- 
2-methyl-4-nitroimidazole (Paul, Kubicki, Jelsch et al, 
2011), 2-methyl-4-nitro-l-phenyl-l//-imidazole-5-carbonitrile 
(Poulain-Paul et al, 2012; Paul, Kubicki, Kubas et al, 2011) and 
for the solid solution of l-(4 / -chlorophenyl)-2-methyl-4-nitro- 
l//-imidazole-5-carbonitrile (97.5%) with 5-bromo-l-(4 / - 
chlorophenyl)-2-methyl-4-nitro-l//-imidazole (2.5%; Poulain 
et al, 2014). 

After high-resolution crystal structure determination and 
multipolar refinement of l-(2 / -aminophenyl)-2-methyl-4- 
nitroimidazole, 1, at 100 K (Paul, Kubicki, Jelsch et al, 2011), 
unexpected high residual-density peaks arranged in a 
'shashlik-like' pattern appeared at high-order residual maps 
(sin 0/k > 0.7 A -1 ) in the planes bisecting the amino groups of 



two symmetry-independent molecules, and a distorted static 
deformation density was observed for one of the nitro groups 
involved in the weaker hydrogen bonds. Thus, third-order 
ANMs were used to model the two fragments of the molecules 
(a split-atom refinement did not succeed). Such a procedure 
resulted in virtually featureless residual electron density maps 
and symmetrical arrangement of the static electron density of 
the N0 2 fragment. 




In the next step powder diffraction experiments at different 
temperatures (20 < T < 300 K) were performed. The data 
collected revealed an isomorphic phase transition (see for 
example Bendeif et al, 2009) as reflected by an abrupt change 
of the b unit-cell parameter around 65 K. Forbidden reflec- 
tions in P2 1 /c did not appear, which suggests that the space 
group was conserved. 

The aim of this paper is an attempt to relate this anhar- 
monic refinement to the isomorphic phase transition by 
analysing several additional X-ray single-crystal diffraction 
experiments performed for 1, including a high-resolution full 
data collection at 10 K on an Agilent Technologies SuperNova 
diffractometer, accurate full data collections at 35 and 70 K 
using the homemade mini-goniometer system implemented on 
an Orange top-loading cryostat on the CRM 2 Bruker AXS 
APEX II diffractometer (Fertey et al, 2007), and temperature- 
dependent unit-cell parameter determination from powder 
diffraction patterns collected on a Panalytical X'Pert Pro 
diffractometer. As careful crystal structure and topological 
analyses of the electron density have already been performed 
for the title compound (Paul, Kubicki, Jelsch et al, 2011; 
Kubicki & Wagner, 2008), they are not repeated in this paper. 

2. Materials and methods 

2.1. Experimental details of X-ray single-crystal diffraction 
measurements 

A yellowish transparent cube-shaped crystal (0.20 x 0.17 x 
0.13 mm) was chosen for data collection at 10 K on an Agilent 
Technologies SuperNova four-circle diffractometer equipped 
with a CCD detector. The temperature was controlled with an 
Oxford Cryosystems cooling device. A total of 2970 frames 
were collected in 35 runs to obtain the high redundancy data 
and 32 additional reference frames were measured to verify 
the stability of the crystal. Diffraction data up to sinO/k = 
1.10 A -1 were collected using the &>-scan method with a 
rotation width of co = 1°. Different exposure times were chosen 
depending on the 20 settings of the detector: 5 s for 20 = 1.25°, 
and 20 s for 20 = -65.45 and 67.95°, with a 55 mm crystal-to- 
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Table 1 

Experimental details for single-crystal measurements at 10, 35 and 70 K - all refinement parameters are given for the multipole model. 



Crystal data 
Chemical formula 
M r 

Temperature (K) 

Wavelength (A) 

Crystal system, space group 

a, b, c (A) 

P(°) 

v(A 3 ) 
z 

D x (g cm" 3 ) 

^000 

Absorption coefficient (mm -1 ) 
Crystal to detector distance (mm) 

Data collection 
Absorption correction 
T T 

1 mm? 1 max 

Crystal size (mm) 
sin 6/X range (A -1 ) 
Limiting indices 

Reflections collected, unique, unique 
with a cut-off 

^int 

Data completeness (%) 
No. of parameters 



C 10 H 10 N 4 O 2 

436.4 

10 

0.71073, graphite-monochromated 
Monoclinic, Fl^ic 

11.0104 (3), 10.0398 (2), 18.6040 (4) 
97.320 (2) 
2039.77 (8) 
8 

1.42 
912 
0.104 
55 



Analytical 
0.983, 0.989 
0.20 x 0.17 x 0.13 
0.07-1.10 

-24 h -> 25, -22 k -+ 22, 
-41 l 41 

305 420, 22 731, 15 217 [/ > 2a(I)] 



0.059 

100 

945 



35 



10.9784 (14), 10.0056 (13), 18.488 (3) 
97.223 (4) 
2014.7 (5) 

1.44 

0.105 
40 



Multi-scan 
0.915, 1.105 
0.12 x 0.14 x 0.15 
0.07-0.90 

-19 -> h -> 16, -15 -> k -> 17, 

-32 I 32 
41 665, 11 032, 9475 [/ > 1.25cx(/)] 

0.065 

96.7 

945 



70 



11.0470 (12), 10.1293 (11), 18.652 (2) 
97.223 (3) 
2070.6 (4) 

1.40 

0.102 
40 



Multi-scan 
0.932, 1.028 
0.12 x 0.14 x 0.15 
0.07-1.20 

-23 h -> 25, -23 -> k -+ 22, 

_44 -> / 44 
121 651, 26 563, 17 738 [/ > 2a(I)] 

0.087 

88.5 

995f 



Refinement 
Weighting scheme 
Goodness of fit on F 2 
Final R(F) indices^ 
Ap max , Ap min (e A" 3 ) 



^=« 2 (\F 0 \ 2 ) 
0.90 v 7 
R 2 = 0.032, wR 2 = 0.028 
0.32 (6), -0.34 (6) 



0.92 

R x = 0.029, wR 2 = 0.028 
0.25 (6), -0.30 (6) 



0.95 

R 1 = 0.029, wR 2 = 0.028 
0.29 (6), -0.27 (6) 



t Different number of refined parameters due to additional ANMs required only at 70 K. t R l = J2\F 0 ~ F c \/J2^o'^ = y^X! [C^o — F c )/aF 0 ~\ 2 /J2 [Fq/ctF^ 2 . 



detector distance. Details of the data collection and crystal- 
lographic statistics are collected in Table 1. 

Another yellowish crystal (0.12 x 0.14 x 0.15 mm) from the 
same crystallization batch was chosen for the 35 and 70 K 
measurements on a CRM 2 Nonius Kappa CCD diffractometer 
equipped with a homemade universal low-temperature mini- 
goniometer, helium top-loading Orange cryostat (Fertey et al, 
2007). A total of 4074 (35 K) and 8240 (70 K) frames were 
collected in 8 (35 K) and 25 runs (70 K). Diffraction data up to 
sin(9/A = 0.90 A" 1 (35 K) and sinO/X = 1.20 A" 1 (70 K) - the 
lower resolution for the 35 K measurement resulted from time 
and liquid helium restrictions - were collected using co- and (p- 
scan methods with 0.25° rotation widths, the x angle being 
fixed at 43.37° and the crystal-to-detector distance at 40 mm. 
Two 20 positions (—30 and —90°) were used to collect all the 
reflections with exposure times of 3 and 60 s for the 35 K data, 
and 10 s and 60 s for the 70 K data. Details of the data 
collections and measurement statistics are given in Table 1. 
Despite some geometrical constraints due to the cryostat 
orientation, the completeness of the data is very close to 
100%, and the internal agreement factors are very good 
compared with typical charge-density quality X-ray data. 

Unit-cell parameter determination, integration of the 
reflection intensities, data reduction and Lorentz-polarization 
corrections were performed using CrysAlis PRO (Agilent 
Technologies, 2013) for the 10 K data, and APEX2 (Bruker, 
2012) for the 35 and 70 K data. An analytical numerical 



absorption correction using a multi-faced crystal model (Clark 
& Reid, 1995) was applied to the 10 K data, while a multi-scan 
absorption correction (Blessing, 1995) was applied to the 35 
and 70 K data. Data sorting, scaling and merging of reflections 
were performed with SORTAV (Blessing, 1997, 1989, 1987) 
for all three datasets. 

As shown in Table 1, all the multipolar models (10, 35 and 
70 K) converge to very good R factors. This shows the possi- 
bility of collecting accurate charge density data using the mini- 
goniometer and cryostat system (Fertey et al., 2007), i.e. 
performing very low-temperature high-resolution accurate 
X-ray data collections with very small helium consumption. 
One of the problems not yet resolved for the mini-goniometer 
data is the precision of the cell parameters (Table 1) possibly 
due to the difficulty in centering the crystal inside the cryostat, 
and to the anisotropy of the data collection needed to avoid 
possible collisions; as shown below, this has some conse- 
quences on the quality of the bond distances and angles. 

2.2. Powder diffraction measurements (PXRD) 

All PXRD measurements were performed using a Panaly- 
tical X'Pert Pro diffractometer equipped with a Cu tube, a 
Ge(lll) incident-beam monochromator (A = 1.5406 A) and an 
X'Celerator detector. Temperature-controlled diffractograms 
were collected on cooling with an Oxford Cryosystems cryo- 
stat (Phenix) from 300 to 125 K (under vacuum, cooling rate 
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6 K min - ; 25 K increments; temperature stabilization: 5 min), 
then from 120 to 15 K (under vacuum, cooling rate 6 K min -1 , 
5 K increments, temperature stabilization: 5 min). Tempera- 
ture-controlled diffractograms were collected on heating from 
20 to 120 K with the same cryostat and then from 125 to 300 K 
under the same conditions. 

Data collection was carried out in the scattering angle range 
0 = 5-50° with a 0.0167° step over 4 h. The program GSAS/ 
EXGUI (Toby, 2001; Larson & Von Dreele, 1994) was used for 
the Le Bail extraction in space group P2 x /c. Owing to the 
complexity of the structure and since powder X-ray diffraction 
(PXRD) is less sensitive than single-crystal measurements, 
single-crystal atomic parameters were used as the structural 
model. Only the cell dimensions, parameters of the pseudo- 
Voigt profile shape function and the zero shift were refined. 

2.3. Structure determination and refinement 

Crystal structures of 1 for the three datasets (10, 35 and 
70 K) were solved using SIR92 (Altomare et al, 1993) and first 
refined with SHELXL (Sheldrick, 2008) applying the inde- 
pendent atom model (I AM), with isotropic and anisotropic 
treatment of H and non-H atoms, respectively. Geometry 
constraints (C Ar -H = 1.083 A; C Me -H = 1.059 A; N-H = 
1.009 A), atomic thermal motion parameters (initial values of 
t/? o = y x Uf q ; y = 1.2 for Ar and NH 2 groups; y = 1.5 for Me 
group) were initially imposed on H atoms to preserve the 
physical meaningfulness of the models. Fig. 1 shows the two 
symmetry-independent molecules of 1 with labelling scheme 
(see Paul, Kubicki, Jelsch et al, 2011, for more details). 

Subsequently the multipolar refinement strategy previously 
presented was applied, with restraints on symmetry and 
chemical equivalency defined as optimal from R free -factor 
calculations (see Paul, Kubicki, Jelsch et al, 2011, and refer- 
ences therein). The main points of the refinement strategy 
were the following: (a) scale factor refined continuously with 
all parameters; (b) anharmonicity parameters refined against 
the high-order data (sin 0/X > 0.7 A -1 ; only deemed necessary 
for the 70 K data); (c) thermal motion and positional para- 
meters for non-H atoms against high-order data alternatively 




Figure 1 

Anisotropic ellipsoid representation of the two symmetry-independent 
molecules of 1 with atom-labelling scheme. Ellipsoids are drawn at the 
50% probability level, H atoms are depicted as capped sticks 
(MERCURY; Macrae et al, 2008). The labels of the second molecule 
are ordered in the same way and marked with an A (e.g. CIA, NLA etc.). 
The strongest interactions are indicated by turquoise dashed lines. 



with H-atom coordinates and distances constrained to stan- 
dard neutron values (Allen & Bruno, 2010); (d) refinement of 
multipolar parameters followed by valence populations 
(constraints imposed on chemically equivalent atoms in a 
similar environment) and then both together; (e) k parameters 
for non-H atoms (constraints imposed on chemically equiva- 
lent atoms in similar environment); (f) points (d) and (e) 
performed until convergence is achieved; (g) anharmonicity 
parameters (only for 70 K data) alternatively with thermal 
motion and positional parameters for all atoms against all data 
(H atoms still constrained); (h) valence and multipole popu- 
lations alternatively with k for non-H atoms and positional 
parameters plus thermal motion; (i) anharmonicity para- 
meters (only for 70 K data); constraints on valence and 
multipole populations together with k, k' coefficients changed 
into restraints at the R free level; (/) k' non-H atoms alter- 
natively with valence and multipole populations; (k) k of H 
atoms; (/) points (h) and (/) alternatively; (m) k' for H atoms; 
(n) point (/) repeated; (n) SHADE estimation of the thermal 
motion of H atoms (Madsen, 2006; Madsen et al, 2013); (o) 
valence and multipole populations alternatively with k H 
atoms, k non-H atoms, coordinates and thermal motion; (p) 
point (j) repeated; (q) point (h) repeated; successive refine- 
ment of k non-H, k H-atoms, k' non-H, k! H atoms; (r) final 
simultaneous refinement of all parameters. 

As mentioned above, only the 70 K data required the third- 
order anharmonic corrections (Kuhs, 1992; S0rensen et al, 
2003) for a correct modelling of three atoms of one N0 2 group 
(N81, 081 and 082) and two amino N atoms (N6 and N6A) in 
order to reduce the typical 'shashlik-like' pattern usually 
found at high-order residual density maps. Stronger interac- 
tions, in which the second nitro group (N8L4, 08L4 and 
082A) is involved, seem to restrict vibrations and therefore an 
harmonic model was sufficient. 



3. Results and discussion 
3.1. Powder diffraction data 

Along with the temperature decrease from 300 to 100 K, a 
linear evolution of the unit-cell volume is observed reaching a 
minimum at ~ 60-65 K (Fig. 2), followed by a slight volume 
increase from 50 to 20 K. The b parameter decreases linearly 
from room temperature (RT) to 60 K with the temperature 
(T) (Ab(T)/b RT = -3 x 10~ 5 T + 0.99) and then increases for 
T < 60 K (Ab(T)/b RT = 2 x 10 5 T + 0.99), in agreement with 
previous findings (Paul, Kubicki, Jelsch et al, 2011). The c 
parameter remains almost constant, as already noted by Paul, 
Kubicki, Jelsch et al (2011). Contrary to the observation of 
Bendeif et al (2009) no hysteresis phenomenon was found or 
if it exists the temperature difference is smaller than 5 K. 

When increasing T in the range 15-100 K, a splitting of the 
100 and 200 reflections appears, which may suggest a second- 
order phase transition (Fig. 3). The phenomenon is more 
pronounced at 100 K, while at 300 K the diffraction peaks are 
practically symmetrical. Such a splitting is not visible on the 
020 reflection due to its very small intensity (Fig. S2). Two 
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Table 2 

Anharmonic nuclear motion parameters greater than 3a for the 100 K and 70 K data. 
100 K 



Cm 


N6 


-0.001724 (81)f 


N6A 


0.000657 (47)t 


N8 




081 


0.000787 (60)t 


082 


0.000494 (66) 


C-222 


N6 


0.000264 (49) 


N6A 




N8 


-0.001359 (69)f 


081 


0.000246 (63) 


082 


0.000379 (95) 


C333 


N6 




N6A 




N8 


-0.000232 (10)t 


081 


-0.000044 (9) 


082 


-0.000034 (8) 


Cll2 


N6 


0.002705 (157)t 


N6A 


-0.000538 (93) 


N8 


-0.000386 (91) 


081 


0.000655 (138) 


082 


-0.000725 (165) 


C\22 


N6 


-0.001220 (128) 


N6A 


0.000432 (92) 


N8 


0.000902 (119) 


081 


0.000411 (141) 


082 


0.000856 (192) 


Cll3 


N6 




N6A 


0.000268 (51) 


N8 


-0.000303 (50) 


081 


0.000324 (71) 


082 




Cl33 


N6 




N6A 




N8 


0.000295 (34) 


081 




082 




C223 


N6 




N6A 




N8 


-0.002084 (85)f 


081 




082 


-0.000378 (101) 


C233 


N6 




N6A 




N8 


-0.001178 (45)t 


081 




082 


0.000386 (44) 


C"l23 


N6 




N6A 




N8 


0.000991 (102) 


081 


-0.000511 (114) 


082 





70 K 



Cm 


N6 


-0.000475 (56) 


N6A 


0.000135 (39) 


N8 


0.000313 (41) 


081 


0.000466 (48) 


082 


0.000519 (52)t 


C222 


N6 




N6A 


0.000422 (34)t 


N8 




081 




082 




C333 


N6 




N6A 


-0.000061 (4)t 


N8 




081 




082 




C112 


N6 


-0.000874 (81)f 


N6A 




N8 


0.000245 (58) 


081 


-0.000604 (76) 


082 


0.000770 (90) 


C122 


N6 


-0.000247 (64) 


N6A 




N8 


0.000321 (51) 


081 


0.000308 (70) 


082 


0.000713 (87) 


C113 


N6 




N6A 




N8 


0.000175 (33) 


081 


0.000348 (42) 


082 


0.000237 (42) 


C133 


N6 




N6A 


0.000095 (17) 


N8 


0.000081 13) 


081 


0.000100 (18) 


082 


0.000064 (17) 


C223 


N6 




N6A 


-0.000627 (36)t 


N8 


0.000113 (23) 


081 


0.000210 (33) 


082 




^233 


N6 




N6A 


0.000361 (18)t 


N8 




081 


-0.000084 (16) 


082 


-0.000136 (17) 


Cl23 


N6 




N6A 


-0.000435 (51) 


N8 


0.000123 (38) 


081 




082 





t 10<T level. 



crystal phases seem to coexist along a large temperature range. 
This can explain the observed disorder at 100 K, which was 
solved using anharmonic atom treatment. 

3.2. Charge-density distribution modelling 

According to our previous findings (Paul, Kubicki, Jelsch et 
al, 2011) for the 100 K data, the largest residual peaks in the 
residual density Fourier maps (sin#/A < 0.9 A -1 ) lie in the 
planes bisecting the H61— N6— H62 moiety, at a distance of ca 
0.5 A from the N atoms (0.37 e A -3 for N6A, and 0.28 e A -3 
for N6 atoms). They disappear at a resolution of sin#/A < 
0.7 A -1 and therefore cannot be interpreted as missing H 
atoms, because they only appear when high-order reflections 



1.005 




-c 
-V 



Figure 2 

Unit-cell parameter variation with temperature decrease from 300 to 
20 K normalized to 300 K. 



are included, while H atoms scatter at very low sinO/k. The 
refined third-order anharmonic parameters are statistically 
not significant, but reduce substantially the residual peak 
heights. 

The first important result of this report is that the 10 and 
35 K data do not need any anharmonic motion modelling 
(ANM) of both amino and nitro groups, whereas ANM 
refinement is still necessary at 70 K as peaks and holes in the 
'shashlik-like' pattern appear close to the N6 atom: +0.42 (6) 
and -0.32 (6) e A -3 [compared with +0.56 (5) and 
-0.27 (5) e A" 3 for the 100 K data (Paul, Kubicki, Jelsch etal, 
2011)]. The lower resolution of the 35 K (0.9 A -1 ) dataset 
compared with the 10 K (1.1 A -1 ) and 70 K (1.2 A -1 ) ones 
does not affect the detectability of the 'shashlik-like' pattern, 
since such a distortion is already observed at 100 K at 0.9 A -1 



300 K 



100 K 



15 K 




— 1 — 

7.8 



— 1 — 

8.0 



1 1 1 1 

16.0 16.2 



8.2 15 
20 (°) 

Figure 3 

Splitting of the 100 and 200 diffraction peaks with temperature increase 
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cut-off (Paul, Kubicki, Jelsch et al, 2011). Moreover, the 
I/cr(I) cut-off at 35 K was reduced to 1.25 compared with 2.0 
for 10 and 70 K in order to improve the data-to-parameter 
ratio. 

Fig. 4 gives residual density maps obtained after harmonic 
(left panel) and anharmonic (right panel) treatment of the 
70 K data. Similar to the 100 K data, residual peaks at 70 K are 
more pronounced for one of the two amino groups (N6A) and 
mostly at higher resolution (1.2 A -1 ). The residual peaks at 
100 K are slightly higher than those observed for the 70 K 
data. Application of the ANMs of third-order significantly 
reduced the positive and negative residual electron density 
peaks and restored the expected valence-density arrangement 
around O atoms in the N0 2 group. 

Comparison of the third-order ANM parameters for the 100 
and 70 K data is given in Table 2 for C i} - k above the 3cr 
criterion. There is a general trend that the significant para- 
meters at 100 K drop considerably at 70 K (e.g. C m for N6 
and N6A atoms). However, surprisingly, some parameters 
seem to be significant only at 70 K (e.g. C 222 , C 223 and C 233 for 
N6A). 

The quality of the four (10, 35, 70 and 100 K) data refine- 
ments is comparable, with insignificant differences between 
the corresponding agreement factors (Table 1): R x = 0.029- 
0.032, wR 2 = 0.025-0.028 and S (goodness-of-fit) = 0.90 (10 K)- 
1.07 (100 K), and Ap max (from +0.25 e A" 3 to +0.32 e A" 3 ), 



p min (from -0.22 e A -3 to -0.34 e A -3 ), which in fact depends 
on the data collection resolution (lowest for 35 K data). 

In conclusion, diffraction experiments at 35 and 10 K did 
not require any special anharmonic treatment, as the 
harmonic approximation is sufficient for all the atoms 
concerned (Fig. 5). It is in line with the isomorphic phase 
transition, which occurs around 60 K. ANH modelling of the 
70 and 100 K data enables modelling of the residual density 
accounting for the disorder which may be due to the co- 
existence of both LT and HT crystal phases existing in this 
temperature range. 

3.3. Electron density model validation via topological 
analysis of the covalent bonds 

In order to compare and validate the model correctness at 
different temperatures (10, 35, 70 and 100 K) the covalent 
bond critical points (CPs) of the aryl ring (that should be 
unchanged and prove consistency of these four data treat- 
ments), together with those of the anharmonic fragments, are 
collected in Table SI. In general, the distance between the two 
involved atoms is ~ 0.01 A longer for the 70 K structure, but 
this lengthening is not significant enough to be reflected in the 
respective distances to the critical points and, as seen below, is 
a result of a less accurate estimation of the cell parameters 
derived from the mini-goniometer data. For the C—C bonds 




Figure 4 

Residual and static electron density maps at 70 K after multipole 
refinement in the planes bisecting NH 2 and N0 2 groups (cutoff 1.2 
and 0.9 A resolution) neglecting (left panel) or including (right 
panel) ANMs; contours set to 0.05 e A -3 , blue dashed lines - 
negative contours, red solid lines - positive contours; 1.2 and 0.9 
values indicate the resolution (A -1 ). 
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Table 3 

Summary of the three strongest nitro group interactions at different temperatures. 
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Cp3 


10 
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1.363 


1.002 


0.076 


1.13 
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0.04 
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35 




2.3489 


1.375 


1.005 


0.069 


1.09 


-0.25 


-0.24 


1.58 


0.03 


23.5 
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70 




2.3812 


1.395 


1.010 


0.065 


1.07 
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0.03 


22.6 
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of the aryl ring the total electron density value differences for 
a given bond are < 0.1 e A -3 , about 2a, while the Laplacian 
values are systematically higher for the 70-100 K data, but 
within the usually accepted estimated error (up to 4.0 e A -5 ). 

For bonds involving the anharmonic atoms the total density 
at CP is on average larger for the datasets, which were 
corrected for anharmonic treatment (maximal change 



Figure 5 

Residual electron density and static deformation maps after harmonic 
modelling of 35 and 10 K data drawn in the planes bisecting both amino 
groups and one nitro group prone to geometrical distortion; contours set 
at 0.05 e A -3 , blue solid lines - positive contours, red dashed lines - 
negative contours, smO/k < 0.9 A -1 . 



0.2 e A 3 for the N8 — 081 bond), as well as the Laplacian 
values ^ 4-6 e A -5 ) for all the 70-100 K bonds, except N8- 
081. 

Contrary to Zhurov et al (2011) the Laplacian maps (Fig. 6) 
of the nitro group calculated within the harmonic approx- 
imation (not shown here) are indistinguishable from those 
correctly modelled, which suggests a lower anharmonicity/ 
disorder in 1. 

For the three critical points characterizing the strongest 
intermolecular interactions where the N0 2 groups are 
involved, the topological data at different temperatures are 
collected in Table 3. All electron density values decrease when 
the temperature increases, while the Laplacian values fluc- 
tuate rather than show a visible trend. Nevertheless, all these 
changes are insignificant at the 3<j level, as expected on the 
basis of constant intermolecular distances [for example, the 
08L4- • -H62A distance equals 2.028 (10) A]. 

In a recent review, Kamihski et al. (2014) investigated 
structural parameters and charge-density properties in a series 
of 100 K high-resolution datasets of a-oxalic acid dehydrate, 
which reveals that electron density and Laplacian values at 
corresponding CPs for this unique crystal structure vary over a 
small range, even at the same temperature. The standard 
deviations for the total electron density and Laplacian for 
covalent bonds and intermolecular bonds vary between 0.03- 
0.06 e A -3 , 1-7 e A -5 and 0.001-0.03 e A -3 , 1-6 e A -5 , 
respectively, which confirms our above conclusion that 
changes in 1 are statistically insignificant. 

3.4. Accuracy of the bond lengths obtained from the mini- 
goniometer data 

As already shown recently (Jarzembska et al, 2013), it is 
difficult to obtain an accurate orientation matrix with the mini- 
goniometer setup, leading to slightly different cell parameters 
compared with those obtained from powder diffraction data, 
which consequently affects the precision of the bond distances. 
Recalculation of the aryl ring C— C bond lengths for the 35- 
100 K data, using the unit-cell parameters obtained from the 
powder diffraction experiment (second row of Table 4), brings 
a much better agreement (Table SI versus Table 5). The 
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Table 4 

Unit-cell parameters of 1 at different temperatures. 
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Powder diffraction 



a (A) 
b(A) 
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b (A) 
c(A) 
/K°) 



11.0104 (3) 
10.0398 (2) 
18.6040 (4) 
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11.0470 (12) 
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97.223 



11.030 (2) 
10.092 (2) 
18.637 (3) 
97.24 (2) 
11.0532 
10.1303 
18.6769 
97.175 



maximal difference in the d n value between 35 and 100 K is 
0.008 A and a clear trend is found: d 12 (35 K) > d 12 (70 K) > d 12 
(100 K). This behaviour has been known for a few decades 
(see for example Busing & Levy, 1964; Scheringer, 1980; 
Destro & Merati, 1995), and results from the higher degree 
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Table 5 

Comparison of the C— C bond distances of the aryl ring for 35, 70 and 
100 K data, using the cell parameters transferred from the powder 
experiment. 



T 






D12 


T 






D12 


( K ) 


Atom 1 


Atom 2 


(A) 


( K ) 


Atom 1 


Atom 2 


(A) 


o c 

JJ 


CI 


cz 




JD 






1 /im 


/U 






1 /inn 


/u 






1 ^os 


1 nn 

1UU 








1 nn 

1UU 








q< 

jj 


CI 


Co 


1 /i 1 q 


Jj 


CIA 


CO/l 


1 /1 1 o 
1.41Z 


/u 






1 a 1 n 


7n 
/u 






1 /i no 


100 






1.409 


100 






1.407 


35 


C2 


C3 


1.401 


35 


C2A 


C3A 


1.396 


70 






1.396 


70 






1.393 


100 






1.395 


100 






1.392 


35 


C3 


C4 


1.404 


35 


C3A 


C4A 


1.406 


70 






1.401 


70 






1.403 


100 






1.400 


100 






1.401 


35 


C4 


C5 


1.400 


35 


C4A 


C5A 


1.396 


70 






1.394 


70 






1.392 


100 






1.392 


100 






1.390 


35 


C5 


C6 


1.418 


35 


C5A 


C6A 


1.415 


70 






1.417 


70 






1.414 


100 






1.415 


100 






1.412 



of precision in determining molecular geometry at lower 
temperatures. 

4. Conclusions 

The aim of this study was to show the link between anhar- 
monicity and isomorphic phase transition of a molecular 
crystal. We have shown that ANM corrections improve the 
charge-density model above the phase transition temperature, 
whereas a simple harmonic model is sufficient below the 
transition temperature. Softening of the anharmonicity is 
therefore connected with the transition mechanism. As shown 
from powder diffraction data the nature of the phase transi- 
tion seems to be second order with a coexistence of both 
phases over a large temperature range (40-50 K). As the 
atomic structures of both phases are extremely similar, a split- 
atom model cannot take into account the disorder observed 
on the residual maps which was accounted for using a third- 
order anharmonic treatment. Such an interpretation however 
needs more experiments on other molecular crystals to be 
considered as a general rule. 



10 K 




Figure 6 

Laplacian of the total electron density maps at 100, 70, 35 and 10 K for 
the two N0 2 groups; anharmonic treatment indicated by ANMs marked; 
logarithmic contours; blue dashed lines - positive contours, red solid lines 
— negative contours. 
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